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Abstract 

The formation of bubble defects of the double stranded DNA is treated ac- 
cording to the Lifshits theory of disordered chains. A molecule of the DNA is 
modelled on a harmonic lattice with nearest neighbour interaction, elastic con- 
stants being randomly distributed. The helicoidal symmetry is accommodated 
through a chiral field at sites of the lattice. The number of sites varies from 
100 to 300, corresponding to DNA segmnets of persistence length. We find the 
spectra of elastic eigen-modes that mimic inter-strand excitations of the duplex. 
The frequency distribution shows peaks and valleys at the high-frequency end 
of the spectra, in accord with the general theory. External excitations may lead 
to a parametric resonance that can generate localized modes of the lattice. In 
real life pumping hyper-sound may generate a resonance similar to that studied 
in this paper, and thus result in excitation of inter-strand modes and possible 
formation of bubbles, in the duplex of the DNA. 



Introduction 



There has been remarkable progress in the studies pQ — [5] of defects of the DNA due 
to breaking hydrogen bonds and formation of local denaturation zones, or bubbles. 
It is generally accepted that the phenomenae are important for understanding basic 
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facts of molecular biology, i.e. denaturation, replication, transcription. A number 
of recent papers have been dedicated to working out a theoretical description of the 
phenomenon, see paper [5] and references therein. Generally, the dynamics of bubble 
formation has been considered either within the stochastic, or nonlinear mechanics, see 
@]. In papers [5], [5], [7] the authors suggest to take into account irregular structure 
of the DNA. In fact, the conformation of the DNA as regards positions of base- 
pairs and values of elastic constants describing its elastic properties, suffers serious 
deviations from the picture of ideal double- helix [5]. The problem is to what extent 
the irregularities of duplex may be involved in the formation of inter-strand defects 
of the DNA. 

The answer may be provided by the theory of dynamic excitations in random media 
worked out, years ago , by I.M.Lifshits and his colleagues, [9] - - [13] (see also [H]), 
p"5]). According to Lifshits' theory defects of media serve centers for the localization of 
elastic excitations, see [5] — [T3] . It should be noted that the introduction of disorder 
results in drastic changes in the shape of spectra, as is illustrated in FigsfT] — [2] 

One can synthesise molecules of the DNA with a prescribed sequence of base-pairs, 
disordered or otherwise, and thus obtain a variety of elastic systems having elastic 
properties required for generating various frequency spectra. Thus, one may expect 
that the DNA may serve a means for experimental study of chaotic elastic dynamics. 

In this respect one may suggest that the formation of bubbles in a molecule of 
the DNA is due to inter-strand motions acquiring an amplitude large enough for the 
hydrogen bonds of a molecule to be broken. Of course, the actual process demands 
taking into account non-linear efects. But the later problem is hard to be solved, 
presently. One may only suggest that generating inter-strand modes of large ampli- 
tude, for example in a resonance regime, is a sign that the bubbles run a chance to 
appear. Therefore, it is worthwhile to look for localzied inter-strand modes in the 
dynamics of disordered chains corresponding to a molecule of the DNA. 

The study of spectra of disordered chains, or lattices, has followed mainly two 
paths. A number of papers employed analytical methods, which generally could 
provide only qualitative description of the spectra, see [14]. It was P.Dean, [16] , who 
developed a powerful computer technic for analysis situs of the frequency distribution 
of disordered lattices. In this paper we use the algorithms worked out by Dean for 
the analysis of inter-strand modes of the DNA of a simple qualitative model which 
mimicks the most essential conformational properties of the DNA. 
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Figure 1: Squared frequency spectra for regular one dimensional lattice . 
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Figure 2: Squared frequency spectra for disordered lattice. Elastic constants selected 
randomly from 1 and 3. 



2 Resonance in random helicoidal lattice. 



It is extremely difficult, if at all possible, directly to employ the theory of disordered 
media indicated above for the needs of the DNA. Therefore, one has to resign to 
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Figure 3: Squared frequency spectra for periodic lattice. Elastic constants are equal, 
masses form sequence ...-1-1,5-1-1-1-1-1,5-1-1-1... 

realistic simplifications of the latter. In building appropriate models we are to take 
into account the following basic features: 

1. a molecule of the DNA consists of two strands; 

2. it has helicoidal symmetry; 

3. its base-pairs generally form a random sequence. 

The double-stranded structure of the DNA provides an opportunity for its base-pairs' 
changing relative positions in space, see FigJU This fact is of primary importance for 
describing the internal, or inter-strand, dynamics of the DNA, in which strands of a 
molecule move from their equilibrium positions. If the amplitudes of the motion are 
small enough the duplex structure of the molecule is preserved, even though deformed. 
The two strands are still joined by hydrogen bonds, their relative positions suffering 
only small changes. But there is a need for a more quantitative description of this 
phenomenon. Aiming at a simple picture of a molecule of the DNA, we choose the one 
dimensional interpretation as a ID lattice whose sites n = 1, 2, . . . , N denote base- 
pairs of the molecule, and describe the change in relative positions of base-pairs with 
vectors Y n . Thus, we obtain a very crude picture of a molecule of the DNA, which 
still preserves the essence of two strands and brings forward the helicoidal dynamics 
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of the system 0. It is worth noting that ID structure of the model is in agreement 
with the conformation of a molecule of the DNA of length less than the persistent one, 
that is 300 — 1000A, depending on ion strength of an ambient solvent. We note that 
a model of the DNA that corresponds to a discreet form of non-linear Klein-Gordon 
equation and take into account irregularities of the molecule, is put forward in paper 
[7], but it does not take into account the helicoidal structure of the DNA. In this 




Figure 4: Site n corresponding to a base-pair. Relative position of the bases denoted 
by vector Y norm . Vector Y m od indicates relative positions of the base-pairs in the 
deformed molecule. Vector Y = Y mo d — Y norm describes the dynamics of the base- 
pair within the framework of the model. 

paper, following [17] . we accommodate the helicoidal structure of the duplex with the 
help of a gauge field describing the helicoidal, or helix, structure d • To that end 
we consider a local frame at each site of the lattice, which provides a reference for 
local strain caused by displacements of adjacent base-pairs. The frames are rotated by 
matrix R due to the transition from cite n to site n+1, and by the inverse matrix 
for the inverse transition from n + 1 to n as is illustrated in FigJS] A deviation from 
the rule results in straining the system and increasing the elastic energy. To provide 
an analytical framework for the above picture we assume that the kinetic energy is 

1 One may recall the familiar story "Dancing Men" by Arthur Conan Doyle, in which crude pictures 
of human beings served a code for the plot. 

2 The concept is quite familiar in gauge field theories, see H.Kleinert, Gauge Fields in Condensed 
Matter, World Scientific, Singapore (1989), and is essentially the same as the angular velocity of 
a top. It dates back to Poisson's time and the theory of motion of rigid bodies, see E.J.Routh, 
Dynamics of a System of Rigid Bodies, Ch. V, Macmillan, London (1892). Similar considerations 
are employed in deriving the Kirchhoff equations for the motion of a rigid body in ideal fluid, see 
H.Lamb, Hydrodynamics, Cambridge University Press, London(1975). 
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Figure 5: Local frames due to the helical symmetry. Coordinate systems x — y and 
x 1 — y', corresponding to sites n and n + 1, transformed into each other by rotation 
R n . The change in values of vectors Y n and Y n+ \ allows for rotation transformation 
R n , according to expression R n +\ — R n ■ Y n . Deviations from the equilibrium position 
coincide with those corresponding to a helicoidal lattice formed by vectors Y n . In 
what follows R n = R. 



related to the dynamics of the vectors Y n and can be cast in the form 



N 



71=1 



-(^) 2 



in which masses m n describes inertial effects accompanying the relative motion of 
bases inside the pairs. An expression for the potential energy should describe both 
the energy related to mutual displacement of the strands, which we cast in the form 



N 



Udev = J2f(Y n f 



n=l 



and the energy caused by the rotation of adjacent base-pairs, which we may write 
down as follows 



U r 



N , 

Et( 4 



+i — R ■ Y n 



n=l 



The total potential energy U is equal to 

U = Urot + U de 
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It is worth noting that Udev equals to zero if Y n — 0, n = 1, 2, . . . , N, as well as 
U ro t = 0, when vectors Y n+ \ and Y n are transformed into each other by rotation R. 
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Figure 6: Localized excitation due to a parametric resonance. All the masses equal 
to 1. Elastic constants are 1 and 1.5 distributed with equal probabilities. Localization 
sites 139 - 148, the excitation's frequency 5.63. 
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Figure 7: Localized excitation with frequency 5.40. Parameters of the lattice are the 
same as in Fig|6] 

Thus, we obtain the Lagrange function of our system 

L = T-U 

where the potential energy U is given by the expressions indicated above. The equa- 
tions of motion read as usual, and, to avoid unnecessary formulae, we shall not write 
them down explicitly. It is important that the equations of motion can be cast in a 
more tractable form. To that end we resort to a simple arithcmctic device by using 
complex notations. The point is that vectors Y n being perpendicular to the axis of the 
molecule we may consider them two-dimensional, Y n = (Yi,!^) and assign complex 
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number z n = Y„ + iY% to vector Y n . Then the rotation given by matrix R introduced 
above corresponds to the multiplication by complex number e l6 ,8 — n/5. Let us 
introduce new variables y n defined by the equation 

y n = e- l& z n , 6 = (n-l)0, n = 1, 2, . . . , N - 1 (1) 
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Figure 8: Squared eigenvalues density G(u> ) histogram for a chain with S — 0.5, e = 0.003 
below normal mode corresponded to eigenvalue 4.006082. 

Using the y n we may cast the Lagrangian of our lattice in the form 

JV r 



n=l 



m n , . . * n k n 2 £n I , 2 

— \Vn ■ y n ) - y 1^™+! ~ y«l ~ y If" I 



where * signifies complex conjugation, and the equations of motion for y ni y n in the 
form 

m.y, = h-xyi-i + k,y l+1 - (k, + - e^i (2) 

and 

miiji* = h^ 1 y*_ 1 + hy* +1 - (k t + - e l y* (3) 
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Figure 9: Squared eigenvalues density G(lo 2 ) histogram for a chain with 8 = 0.5, e = 0.06 
below normal mode corresponded to eigenvalue 4.012972. 



Thus we split the equations of motion for the hclicoidal lattice into two complex 
equations (HI),© for a simple harmonic lattice, the chiral field given by matrix R 
being accommodated by substitution ([T]). Now we may use the method worked out 
by Dean, [16j , for harmonic chains. The fact that the chain variables, y n , are complex 
numbers is not essential. Details of Dean's algorithm see in the review article |16j . 

At this point it is necessary to indicate a feature important for the problem at 
hand. It is generally accepted [TS], [IB], 0, that eigenmodes become localized "even 
for a very small chaotic perturbation of the interaction potential" , [5], p. 39. The state- 
ment does not make any claims as to the size of localization, but our calculations 
indicate that a magnitude of disorder given by values of elastic constants is essential. 
We have made a series of calculations in which the elastic constants take on values 
1 or 1 — e with equal probability, and found that depending on S, the spectra of 
eigenmodes and the sizes of localization acquire the following shapes: 



1. S = 0.001 - distribution function does not change appreciably; no localization 
on the scale of the lattice; 
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Figure 10: Squared eigenvalues density G(u> 2 ) histogram for a chain with 8 = 0.5, e = 0.009 
below normal mode corresponded to eigenvalue 4.020306. 



2. S = 0.006 - distribution function does not change appreciably; no localization 
on the scale of the lattice; 

3. S = 0.009 - distribution function looks the same but there is localization on the 
scale of the lattice, see 

4. S = 0.03 - distribution deformed; localization 30 - 40 sites involved; 

5. S = 0.05 - distribution deformed; localization 30 - 40 sites involved; 

6. S = 0.06 - distribution deformed; localization within a small region. 



Therefore, we need to take into account values of random disturbance of elastic con- 
stants, if we are looking for localizations of a size less than that of the sample. In a 
situation pertinent to a molecule of the DNA, it is reasonable to assume that elastic 
constants may have values of relative proportions 1 and 1.5, in accord with the num- 
ber of hydrogen bonds, 2 and 3, for the base-pairs A-T and G-C. It enables to obtain 
fairly pronounced localizations confined to small regions of a lattice. 
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Figure 11: Squared eigenvalues density G(u> 2 ) histogram for a chain with S = 0.5, e = 0.05 
below normal mode corresponded to eigenvalue 4.124748. 

The data and figures indicated above suggest that there is a large amount of 
situations in which localized modes could be expected for duplexes of the DNA. One 
may employ parametric excitations to generate them. Taking into account attenuating 
and damping effects we consider the dynamics described by the following equation 

miiji = ki-iyi-i + hy i+ i - (hi + h-i)yi - eiViV - iy% + F t (4) 

and choose external excitation Fi (t) given by 

Ft = 2A jft sin (2ut) (5) 

It is important that there is a parametric resonance for frequency v which is the 
double of that of a lattice eigenmode. The amplitude A of the excitation is to be large 
enough to counterbalance dissipative effects due to the term 72/j. Since masses of 
the lattice sites, as well as dissipative and elastic constants, 74, and fcj, are randomly 
distributed. It is difficult to write down an analogue of the Rayleigh condition, [15] , 
for the resonance 

(w 2 - i/ 2 ) 2 = 4 2 - (I) 2 v 2 
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Figure 12: Squared eigenvalues density G(uj 2 ) histogram for a chain with S = 0.5, e = 0.06 
below normal mode corresponded to eigenvalue 4.168756. 



If all the masses, as well the dissipative and elastic constants, are equal, the above 
equation is valid for an eigenmode with frequency u. At any rate, we may conclude 
that there is a threshold for amplitude A, below which one cannot obtain a resonance. 

One may put forward the following arguments in favour of our approach. The 
interactions within base-pairs being estimated, [19] , less 2fceT, thermal excitations 
of energy ksT are generally accepted as a cause of DNA breathing, i.e. opening and 
re-closing of hydrogen bonds, [T] , [5] . It is important that they admit an interpretation 
using wave theory. L.I.Mandelstam, [52], indicated that waves of density introduced 
by A.Einstein in the theory of Brownian motion, [23 , are elastic waves, so that one 
can consider density fluctuations as a result of interference of the elastic waves, [213] . 
[21] , or hyper-sound. In assessing the part played by hyper-sound in the inter-strand 
dynamics it is necessary to compare characteristic sizes. The diameter of a molecule 
of the DNA is approximately equal to 20A, the length of its segment comprising 
10 base-pairs is 34A, whereas the wavelength of hyper-sound in the range of a few 
tens of GHz is of order of several tens A, that is all characteristic sizes involved 
are of order 10 -6 -j- 10 cm. The sizes are propitious for generating inter-strand 
motions. Indeed, imagine a hyper-sound wave of wave-length approximately 40A, 
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Figure 13: Picture on the left hand side illustrates the regular configuration. The 
central and the right ones show localized excitations, frequency 5.63, localization at 
139-148 base-pairs. 



that is frequency oc 10 12 Hz. The chances are that its crest may be located at one 
strand of the molecule whereas its bottom at the other. Therefore, there will emerge 
a torque and a strain, periodic in time, that pull apart the opposite strands. If 
amplitudes of parts of the molecule arc small enough it is only natural to suggest that 
the motion is to be described by terms linear in fields Y n . Thus we arrive at the form 
|5]) for the external excitation^]. 



3 CONCLUSIONS 



According to the arguments given above inter-strand defects, or bubbles, of a molecule 
of the DNA could be due to an interplay between random sequences of base-pairs and 
hyper-sound. The latter is an agent that generates density fluctuations, [21j that are 
assumed to break hydrogen bonds.. Thus, we may suggest that pumping of hyper- 
sound could make for bubble defects, through maintaining inter-strand modes due 
to the parametric resonance effected by hyper-sound pumping. Whether this process 

3 Since we do not aim at studying the denaturation and other non-linear phenomena, we refrain 
from introducing clever non-linear terms in the potential energy of a molecule of the DNA, and 
confine ourselves to harmonic approximation. 
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really takes place depends on magnitude of dissipation accompanying the inter-strand 
dynamics. Presently, it is often claimed that the attenuation is so strong that the 
inter-strand modes should be overdamped. The arguments are based on estimating 
relaxation due to viscosity due to surrounding liquid medium, 24J . But in estimating 
the dynamics of liquid inside GHz-region one must take into account that there the 
Navier-Stokes equations cease to be valid. [2"T]. In this region one has to turn to the 
theory of GHz-hydrodynamics worked out by Mandelstam and Leontovic, [20], [21] . 
So far there have been no estimates of the attenuation made within the framework of 
their theory. In this respect, it is worth noting that the Navier-Stokes theory gives to- 
tally inadequate results for the treatment of the Mandelstam-Brilloiun light scattering 
due to hyper-sound. In fact, there are different opinions as to the importance of the 
disspative effects: M.E.Davis and L.L. VanZandt, [25], following Maxwell's approach 
to hydrodynamics, came to the conclusion that the attenuation of inter-strand modes 
does not overdamp them; in paper [26] the authors, on different grounds, come to 
the same conclusion. Thus, the question of the attenuation of inter-strands is still 
open. Nonetheless, considerable efforts should be taken to circumvent undesirable ef- 
fects of attenuation in experimental research. The important advaces in hyper-sound 
acoustics, in particular constructing transducers, [28] , enables us to be reasonably 
optimistic as to the experimental feasibility of hyper-sound investigation of the DNA. 
One may suggest that the attenuation of hyper-sound could be a means for detecting 
structural changes in the DNA. It is worth noting that to that end, using orieta- 
tionally organized structures of the DNA, for example its liquid crystalline phases, 
[27] . would be useful. Atomic force microscopy could be an important instrument as 
well for studying the formation of the bubbles, even more so as it permits to employ 
molecules of the DNA on solid substrate, and thus to escape at least difficulties due to 
the disspative effects as regards molecules of the DNA in solution. For this argument 
we indebtted to K.V.Shaitan. 

Finally, we would like to draw attention to the fact that a more refined but still 
practical model of the DNA, for example that worked out in paper [29], could make 
for better understanding chaotic dynamics of inter-strand modes, their bearing on 
physics of the DNA, and possible biological implications. 

We are thankful to K.V.Shaitan for the useful discussion. 

The paper had the financial support of the Russian Fund for Fundamental Re- 
search (RFFI) # 09-02-00551-a. 
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